Prospects for the cavity-assisted laser cooling of molecules 
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Cooling of molecules via free-space dissipative scattering of photons is thought not to be prac- 
ticable due to the inherently large number of Raman loss channels available to molecules and the 
prohibitive expense of building multiple repumping laser systems. The use of an optical cavity to 
enhance coherent Rayleigh scattering into a decaying cavity mode has been suggested as a poten- 
tial method to mitigate Raman loss, thereby enabling the laser cooling of molecules to ultracold 
temperatures. We discuss the possibility of cavity-assisted laser cooling particles without closed 
transitions, identify conditions necessary to achieve efficient cooling, and suggest solutions given 
experimental constraints. Specifically, it is shown that cooperativities much greater than unity are 
required for cooling without loss, and that this could be achieved via the superradiant scattering 
associated with intracavity self-localization of the molecules. Particular emphasis is given to the 
polar hydroxyl radical (OH), cold samples of which are readily obtained from Stark deceleration. 

PACS numbers: 32.80.Lg, 42.50.Vk, 39.10.+j, 32.80.Pj 



I. INTRODUCTION 

The experimental realization of large samples of ul- 
tracold, ground state dipolar molecules would be a ma- 
jor breakthrough for research in fields as diverse as ul- 
tracold collisions and chemistry to quantum information 
processing and the study of novel correlated states of 
matter [IJ. In particular, the anisotropic dipole-dipole 
interaction becomes non-negligible for polar clouds below 
approximately 100 /xK. Exotic states of dipolar matter, 
such as field-linked states and dipolar crystals, may be 
observable in this regime [5]. If the rich field of ultra- 
cold alkali Feshbach physics is any measure, then ultra- 
cold molecular collisions and chemistry in the presence 
of the dipolar interaction and external electrical or mag- 
netic fields promises to be fascinating [HHJin]. Moreover, 
the precision motional control attainable only at ultralow 
temperatures is crucial for constructing the architectures 
necessary to realize quantum logic gates or spin lattice 
simulations using the dipolar interactions [Hl[7]. 

While many techniques for ultracold, ground state po- 
lar molecule production show promise, none so far have 
simultaneously yielded the low temperatures and high 
densities required to pursue these goals. Photoassocia- 
tion of ultracold atoms and subsequent optical pumping 
to the molecular ground state [H [9] has achieved lower 
sample temperatures (~ 100 ^K) than techniques such 
as buffer gas cooling (~ 400 mK) |10j. The Stark decel- 
erator provides a nice compromise between density and 
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temperature 112] . Electric and magnetic trapping of 
samples as cold as 10 mK at densities approaching 10^ 
cm"'^ have been demonstrated |13[ 114] . However, new 
cooling techniques are required if we hope to push well 
below the 1 mK regime. 

Unlike atoms, molecules typically have a large number 
of channels into which a given excited state can decay. 
This makes the efficient free-space laser cooling [T5] of 
molecules challenging due to practical limits on the num- 
ber of lasers one can build and operate to achieve ground 
state repumping after each photon scattering event [TB] 
(for a practical scheme involving only several repumpers, 
see [17]). Evaporative and sympathetic cooling tech- 
niques are quite promising, but require an initial den- 
sity higher than what is currently available and is sen- 
sitive to the particular molecular species' collision cross- 
section, which is generally unknown [18]. Cavity-assisted 
laser cooling [THl HO] is a promising — though not fully 
understood — solution in that it provides dissipative cool- 
ing largely independent of the details of the molecular 
structure when the coupling between the molecules and 
the cavity is strong. That the coupling needs to be strong 
is of prime importance for achieving the practical cooling 
of molecules. Establishing this fact as well as exploring 
how to accomplish strong coupling between the cavity 
and molecular system is the main subject of this Article. 

Cavity-assisted laser cooling of molecules brings to- 
gether many otherwise disparate disciplines: cavity quan- 
tum electrodynamics (QED), molecular physics, and 
laser cooling and trapping. In this Article, we describe 
the physics of cavity-cooling species that do not possess 
closed transitions, and in doing so, we review some el- 
ements of cavity QED and atom-photon interactions so 
as to build a common language amenable to researchers 
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FIG. 1: a) Fabry-Perot cavity with resonance frequency ojc 
driven by a laser of frequency ujd and drive strength Qd, de- 
fined in Eq. [4] In steady-state, there are ric — {a'' a) intracav- 
ity photons which escape the cavity via the mirrors at rate 2k. 
b) Three-level representation of the A'^ intracavity molecules. 
A classical laser field of Rabi frequency, ilp, and frequency, 
LOp, couples the ground state \a) to the electronically excited 
state |e), whose frequency difference is uia- Excitations de- 
cay at the total rate 7. Assuming |e) is unsaturated, pop- 
ulation is elastically (Rayleigh) scattered back to \a) at the 
rate yny. Population is lost to the myriad molecular states, 
represented collectively by |6), at the Raman scattering rate 
7_R„. The difference between uia and the transition frequency, 
uJb, between \b) and |e) is much larger than Apt,, c) The 
molecules-cavity system may be pumped by the cavity drive 
laser, a transverse pump beam (shown), or both. The trans- 
verse pump beam of frequency cup and Rabi frequency, Qp, is 
typically red-detuned from both the cavity and the molecular 
resonance {ujp < ujc <^ uJa)- The A*' molecules may be trapped 
or transiently passing through the cavity mode. 



from all three fields. We begin by reviewing the cav- 
ity cooling of a two-level system. This is followed by 
a discussion of molecules in cavities. Whereas strong 
coupling is not necessary for cooling two- level systems, 
it is shown that it is crucial for cooling molecular sam- 
ples. We conclude by analyzing methods for increasing 
the efficiency of cooling molecules given the constraints 
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FIG. 2: The transverse pump laser (red) scatters coherently 
from the far-detuned atomic medium (green). This laser does 
not excite much population, but its energy distribution can be 
Doppler broadened due to the velocity spread of the molecu- 
lar cloud. The blue-detuned and overlapping cavity resonance 
(dashed) modifies the scattered spectrum to favor higher en- 
ergy photons, thus extracting energy from the atomic motion. 
Figure adapted from Ref. [2T) . 



posed by experimentally realizable cavity geometry, cav- 
ity quality, and sample densities and temperatures. For 
a concrete discussion, we focus our attention on the OH 
radical, a highly polarizable molecule well suited for ex- 
ploring dipolar physics [2, 22^, cold coUisions [5, 23^, pre- 
cision measurements of fundamental constants [Mj [5S], 
and quantum information processing [6l I25j . 

The cavity cooling system is illustrated in Fig. [T] in 
which a standing wave cavity is combined with a three- 
level system. While cavity cooling was originally dis- 
cussed in the context of cooling via driving the cavity 
mode |191 126] , we will focus on pumping via transverse 
beams, which provides strong cooling in all three dimen- 
sions [571 121] • Figure ^c) depicts the configuration for 
3-D transverse pumping, though the two pump beams 
in the plane perpendicular to the page are not shown. 
Throughout this Article, unless otherwise stated, discus- 
sion is confined to the transverse pumping scenario. 

From a semiclassical viewpoint, the molecules act as 
a nonlinear element that coherently transfer excitation 
from the red-detuned transverse pump field to higher- 
frequency (bluer) intracavity field. Because the blue- 
detuned field leaks out of the cavity stochastically, en- 
ergy proportional to the fields' frequency difference is 
extracted from the system at the rate 2k. Consequently, 
the molecules' motion is cooled, as this is the only mech- 
anism from which to extract this energy difference. The 
energetics of this cooling mechanism are shown in Fig. |2] 

The primary purpose of this Article is to demonstrate 
that to make this molecular cavity cooling mechanism 
efficient — and therefore practical — one more condition 
must be satisfied. Namely, the molecules must coher- 
ently scatter photons into the cavity mode at a rate much 
faster than the scattering into free-space. This is impor- 
tant not because the free-space scattering leads to ad- 
ditional recoil heating, but because free-space scatter- 
ing allows the molecules to relax back to the multitude 
of metastable molecular states — spin-orbit, rotational. 
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and vibrational — via inelastic, Raman scattering. This 
prematurely quenches the cooling process because these 
states are effectively dark to the very far off-resonant 
transverse pump laser and will not be repumped. We 
address the conditions that must be satisfied to ensure 
that this Raman loss rate is less than the cooling rate, 
and discuss experimental constraints to these conditions. 



II. SINGLE PARTICLE MASTER EQUATION 
ANALYSIS 



The previous section described the general idea of 
cavity-assisted laser cooling, but greater insight into the 
problem is required to evaluate its effectiveness for cool- 
ing molecules with realistic internal energy level distri- 
butions. Solutions to the master equation for the joint 
molecule-cavity system allow the identification of opti- 
mal experimental conditions for achieving cavity-assisted 
laser cooling. However, to gain intuition, we first study 
the case of a single, motionless, two-level atom cou- 
pled to only one cavity mode. In the following sections, 
we use the term "atom" when discussing general cavity 
QED and cavity cooling, and reserve "molecule" to when 
we discuss phenomena or conditions peculiar to cooling 
multi-level systems. 

The derivation of the interaction-picture Hamiltonian 
begins by applying the unitary transformation U = 
exp ^—iujpt{a^ a + (t+ct-)] to the joint atom-cavity system 
Hamiltonian, which assumes that the drive and pump fre- 
quencies are approximately equal {ud ~ ujp). The electric 
dipole and rotating wave approximations are then applied 
to arrive at the Hamiltonian describing the coherent dy- 
namics in the presence of both the drive and the pump 
fields: 

+np[a-+a+]/2 + n4a + a^]/2. (1) 

In this equation, (t_ is the atomic lowering operator, and 
a is the cavity field annihilation operator. The first two 
terms are the bare atom and cavity energies, with Apa = 
— and Ape = Up — Uc, while the atomic pump and 
cavity drive terms are the fourth and fifth, respectively. 
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The atom-cavity detuning is Aca = A. 

The third term represents the atom-cavity interaction: 
excitation is coherently exchanged at the rate g, which 
depends on the atom's position in the intracavity mode 
structure and is inversely proportional to the square-root 
of the cavity mode volume. Specifically, g = gQ^{f) with 



ng = ^i-£ = i.{f)t,.L^^^{f)n.r^^ (2) 

where ^^{f) < 1 accounts for the atom's position, pi is the 
transition dipole moment, and £ is the electric field. We 
have explicitly included Plank's constant in the above 
equation for clarity, but will set h = \ in all subsequent 



equations, including Eq. [T] The mode volume is Vm ~ 
ttwqL/4, where wq and L are the cavity waist and length, 
respectively. In the limit that the atomic linewidth is 
dominated solely by radiative processes, 7 — 2"f_i_, where 
7_L is the decay rate of the atomic dipole. In a sense, 
2g is the Rabi flopping rate of the atom stimulated by 
the vacuum field of the cavity. The purpose of making 
small cavities in single-atom cavity QED research [29 , 30] 
is to enhance the electric field associated with a single 
photon so that this zero-point field fluctuation has a non- 
negligible effect on the atom. 

When one has an atom in a transversely pumped cav- 
ity, the situation is not unlike that of a free space atom 
interrogated by two lasers with Rabi rates Op and 2gy/n^, 
where ric = {a^a)o is the number of intracavity photons 
in the absence of an atom. The presence of the atom in 
the cavity modifies the number of intracavity photons, 
and this back reaction of the atom onto the intracavity 
field produces a self-consistent field intensity that is not 
necessarily equal to (a^a)o. 

The master equation for the density matrix, p, of the 
joint state of the atom and cavity with a transverse pump 
field is the following: 

p — —i H,p + -f±{2(T^p<7+ — <7+(T-p — pa+(j-) 

+K{2apa^ — ap ~ pa^ a). (3) 
The cavity energy decay rate is 2k, and ric is equal to 



A2 ■ 



(4) 



All optical fields are far detuned from the atomic reso- 
nance in the cavity cooling scenario, and typically Ap^, <^ 

{Apa , Aca] . 

Steady-state analytic solutions to the master equation 
may be obtained in the limit that coherences between 
the atom and cavity vanish on time scales much faster 
than any other in the problem. The resulting semi- 
classical equations do not properly account for quantum 
fiuctuations of joint atom-cavity excitations. Neverthe- 
less, it may be shown that it reliably models the single 
atom-cavity behavior in the limit of: 1) weak driving 
wherein the atom is unsaturated — a condition we will 
later adopt; 2) very strong driving where the atom is 
completely saturated; and 3) when there are several in- 
tracavity atoms. Cases 1 & 2 may be understood from 
the spectrum of joint eigenstates for the atom-cavity sys- 
tem shown in Fig. [3] Weak driving means that the first 
excited "rung" on the ladder is minimally populated, 
while in strong driving, many rungs are populated and 
the excitation spacing becomes approximately equal, i.e., 
y/n^ ~ -^/tIc + 1. In both cases, the dynamics may be 
treated semiclassically. 



We will discuss case 3 separately in Section |VIIH but 
essentially, the uncorrelated nature of the atoms aver- 
ages away the quantum fiuctuations, even in the inter- 
mediate regime between weak and strong driving. The 



4 



Tiuj 



2% 



flU) 



|a,2) + |e, 1) 
|a,2)-|e,l) 

|a,l) + |e,0) 
|a, 1) — |e, 0) 

|a,0) 



FIG. 3: Jaynes-Cummings ladder representing the spec- 
trum of joint atom-cavity eigenstates for the case when the 
atom and cavity are on resonance (wc = u>a)- The kets 
|atom, cavity) denote the joint states. The ground state, 
ja, 0) = \g), has zero excitation while the first band of ex- 
cited states are symmetric and antisymmetric superpositions 
of one excitation in either the atomic |e,0) or cavity |a, 1) 
modes. These so-called dressed states [JT] are coupled by 
the "single-photon" Rabi rate, 2g. The next excited band is 
similar, but with a ^JnZ enhancement of the coupling |29) . 



multiatom case exhibits optical bistability in accordance 
to the solutions admitted by the nonlinear semiclassical 
(optical-Bloch) equations. In contrast, a bilinear (mas- 
ter) equation describes the coupling of a single atom to 
the cavity mode, and thus, does not exhibit nonlinear 
effects in steady-state even though the atom-cavity cou- 
pling may be strong. This is due to the dominance of 
quantum fluctuations in the single atom case that, in a 
sense, "washes-out" the nonlinear effects [32]. We men- 
tion this to highlight the fact that in the far-detuned, 
many particle regime of the cavity-cooling scheme, semi- 
classical approximations are appropriate and will be used 
throughout this Article. 

To obtain the steady-state semiclassical solutions, we 
make the approximation that all mixed operator expecta- 
tions are factorable, e.g., {a)^-) ~ {a^){a-). This allows 
us to separate Eq. [3] into two Hamiltonians: one for the 
atom, Ha, in which the cavity operator is converted to 
a C-number a a; and one for the cavity. He, with 
(7_ C; where a is the field amplitude and C is the 
atomic dipole: 



Ha ~ -Apa(T+a-_ + g[a*a-_ + aa^] + ripl&j^ + &_]/2 
= -Apaa+a-+n'p[a+ + a^]/2 (5) 

He = -Apca'a + g[(a^ + C a] +nd[a + d'^]/2 

(6) 



The second lines use the effective atom drive fi^ = 2ga + 
rip and cavity drive fl'^ = 2g^+fld for clarity. The master 



equation may be similarly separated: 

Pa = -i Ha, Pa + 1 ^{'^CF ^ PaO + - U +Cr ^ Pa - Pa^f +Cr -) 

Pc = -i HcPc + K{2apca) - a)apc- PcO^a). (7) 

Setting Pa = = Pc, we obtain for the steady-state 
cavity field amplitude: 



Tr[pca] 



A 



pc 



and for the atomic coherence: 

%/2{Apa - ill) 



C = Tr[pa(T_] 



\^'p\y2 + ll + Ala 



(8) 



(9) 



from which we obtain the atomic excited state popula- 
tion: 



\n'p\^/2 + il + Ala' 



(10) 



To solve the system of Eqs. [8] and |9] in closed- form, one 
needs to assume that the atomic excited state is not sat- 
urated, i.e., 7^ + ^ Q,'p . Fortunately, this is exactly 
the same weak driving condition we must satisfy to pre- 
vent the atom or molecule from incoherently scattering. 

The saturation parameter, s, characterizes the bound- 
ary between incoherent and coherent scattering: 



nl/2 



(11) 



defined here in the absence of a cavity field. When be- 
low s ~ 0.01, the frequency of coherently scattered pho- 
tons matches the incoming, and they are emitted in a 
dipole pattern with respect to the drive beam. In con- 
trast, a saturated atom scatters photons incoherently into 
Ait with an atomic linewidth broadened spectrum. The 
rate of scattered photons is suppressed in the unsaturated 
regime, as may be seen from the following relation: 



In the limit that s < 1 



7 s 
2{\ + s) ~ 

%/2 



Apa + nj_ 
and the cavity field amplitude becomes: 

gnp/2 + nd{Apa + ij±)/2 



a •■ 



{Apa + iji_)iAp 



IK) 



(12) 



(13) 



(14) 



III. THE ROLE OF COOPERATIVITY 

For the discussion in this section, we will set the cav- 
ity drive to zero, fid = 0, and discuss cavity cooling via 
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transverse pumping solely, as sketched in Fig. [TJc). We 
will address fid ^ in a, subsequent section. Energy 
enters the system only through the pump beam, whose 
intensity is proportional to flp. Energy leaves the system 
either via the cavity mirrors at a rate proportional to 
the intracavity field, Tc — 2K|ap, or by relaxation of the 
atomic excited state to modes other than the cavity, at 
rate Ta- We assume that the solid angle subtended by the 
cavity mirrors is much less than 47r. Indeed, for all cavi- 
ties of interest, the fractional solid angle is ^ 10~^, and 
Ta is effectively unmodified by the cavity's presence. In 
light of the discussion in the Introduction, Tc is propor- 
tional to the cooling rate and therefore desirable to en- 
hance, whereas the photons scattered at the rate Ta lead 
to recoil heating and possibly to inelastic transitions — as 
in the case of molecular systems — and is therefore to be 
minimized. 

To ensure low atomic saturation, we now apply the 
dispersive weak driving condition. 



to arrive at the following rate expressions: 



4A2 A2 - 

pa pc 



Ta ~ 



4A2/ 



(15) 

(16) 
(17) 



In the previous expressions, we have made the additional 
assumption that Ap^ 3> Ape, which is true for optimal 
cooling, which in turn requires Ape ~ —k. In light of the 
previous statement, Eq. [16] equals 



2k 



(18) 



and the ratio of scattered photons that contribute to cool- 
ing to those contributing to heating (and Raman loss in 
the case of molecules) is: 



C 



Ta 
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2k7_l 



(19) 



The label "C" is chosen to highlight that this enhance- 
ment factor is none other than the single-atom cooper- 
ativity parameter well-known from cavity QED |29j . In 
essence, the condition C 1 signifies that the coher- 
ent atom-cavity dynamics dominate over dissipative. In 
other words, the cavity field is stimulating the transfer 
of photons from the pump field to the decaying cavity 
mode at a rate faster than free-space scattering. The 
cooperativity may be better understood by noting that 
the decay rate of the coupled atom-cavity system, in this 
hmit, is Tac « 7(1 + 20)^1/ Aja EH]- The extra factor 
2C comes from the increased rate at which the pump 
field is depleted. Important to note is the independence 
of the cavity-to-free space scattering ratio, C, on detun- 
ing. Contrary to expectations, one cannot escape the 



free-space scattering that produces Raman loss by de- 
tuning far from resonance. Detuning only leads to slower 
cooling, since both rates depend on the saturation pa- 
rameter in the same fashion: once incoherent scattering 
is quenched, there is no further advantage to increased 
detuning or decreasing pump intensity. 



IV. CAVITY COOLING IN THE 
PERTURBATIVE LIMIT 

When C < 1 (weak-coupling limit, assuming g < 
['*5 7±])i the cooling process may be described pertur- 
batively as coherent Rayleigh scattering, as was done by 
Vuletic et al. [151 HZ]. A remark on language: "coherent" 
refers to scattering in the s < 1 regime, wherein a definite 
relationship exists between the scattered field and the os- 
cillating atomic dipole. "Rayleigh" refers to the fact that 
any light scattered into the cavity mode is of a frequency 
such that the atom relaxes to the ground state of origin, 
and is thus an elastic process. This occurs when Ape 
is much smaller than any frequency difference, cOab, be- 
tween I a) and metastable ground states, |&). The phrase 
"coherent Rayleigh" means that the Rayleigh scattering 
is in a dipole pattern and is of the same frequency and 
spectral bandwidth as the pump field. 

For consistency, we remark that while the notion of 
coherence in Raman scattering — which is inherently an 
inelastic process — is not commonly defined, we take the 
phrase "coherent Raman scattering" to mean that the 
character of the scattered field retains a definite relation- 
ship to the pump field, which is the case in the s 1 
regime. While there is a frequency offset between the 
pump and scattered field, there is a fixed phase rela- 
tionship between them, and the field is scattered in a 
dipole pattern. As will be discussed extensively in Sec- 
tion |VT] this coherent Raman scattering, which quenches 
the molecular cooling process, cannot be overcome by 
making s small. The only remedy is to enhance the 
molecule-cavity coupling. 

Since 5 < k in this limit, "scattering" is the term of 
choice because any photon emitted into the cavity es- 
capes via the mirrors before being reabsorbed by the 
atom. Thus, the cavity field is incidental: the atom 
scatters photons from the incoming pump beam to the 
two beams emanating from the cavity mirrors and the 
cavity itself serves only to provide a concentrated den- 
sity of states, which modifies the frequency spectrum of 
scattered photons. The language of the Purcell effect — 
enhanced emission into the solid angle subtended by the 
cavity mirrors — is apt in this situation and Vuletic et al. 
makes use of it to describe these dynamics. As shown 
in Appendix |Aj the Purcell factor is exactly equivalent 
to the cooperativity, but expressed in a more experimen- 
tally recognizable form. We will return to this formalism 
when discussing multimode cavities in Section |VII[ 

The cooling force and rate can be readily derived in 
the perturbative regime, as was done by Vuletic et al. 
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for the case of cavity pumping (fi^ ^ 0, flp = 0) [T^] and 
transverse pumping (fi^ = 0, fip ^ 0) [22] • This is done 
without invoking the approximation kv < k, where k is 
the optical transition wavenumber and v is the atomic ve- 
locity. Wc do not revisit these derivations, but merely de- 
scribe the conditions for obtaining optimal cooling rates 
and forces for the transverse pumping case. 

In the coherent Rayleigh scattering process, the energy 
extracted is equal to the difference in energy between the 
cavity resonance and the frequency of the pump beam. 
The rate of cooling and damping force, however, depend 
on the spectral overlap between the scattered pump field 
and the cavity resonance. Figure |4] shows the coordinate 
system for cooling in the 2D plane spanned by one set 
of transverse pump beams and the cavity axis. Cooling 
occurs via a two-photon scattering process in this pertur- 
bative limit, with a lower energy photon absorbed along 
Hkx and a higher energy one emitted along itk^. The 
cooling force is thus along the diagonals spanned by k^^ 
and kz. 

The scattering rate into the cavity from a single trans- 
verse beam(e.g., -t-k^,) may be found by incorporating the 
Doppler shift k • v into the Ape detuning parameter of 
Eq. 16 Note that we assume \Apa\ ^ \kv\. Making the 



substitution 



|Apa| > , , 

• V, we have: 



rc(v) = 2k 



4A2 

pa 



(20) 



In the above expression, we have ignored the recoil- 
induced detuning, ujrec — hk'^ /2m, since this is 100 times 
smaller than the Doppler shift for the ^10 m/s OH sam- 
ples we consider here. The cooling power from this beam 



Pcooling OC HkvT , 



(21) 



and the cooling rate for a velocity, v, that is much larger 
than the recoil velocity hk/m, is: 




FIG. 4: Illustration of the cavity cooling mechanism present 
in the perturbative (weak coupling) limit. This mechanism 
can be referred to as cavity Doppler cooling [T^ in analogy to 
free-space Doppler cooling. A scatterer preferentially absorbs 
red-detuned photons from the +kx beam. Its forward veloc- 
ity in blue shifts emitted photons into resonance with 
the blue-detuned cavity, enhancing the scattering rate in this 
direction. As a consequence of this two-photon process, the 
scatterer receives a momentum kick in a direction opposite 
to its combined and velocities and is therefore cooled. 
A similar process works for cooling in y with the addition of 
another pump laser pair. Figure adapted from Ref. |27| . 



cavities and light molecules such as OH, the recoil limit. 
Tree = hojree/kB, Can be larger than the cavity cooling 
limit, Tf. The final temperature will therefore be limited 
by recoil. For two- level atoms, we see that one does not 
need to be in the strong coupling regime to achieve cool- 
ing. Cooperativity less than unity only slows the cooling 
rate and prevents the system from achieving the lowest 
possible temperature, set by Hk. This no longer holds 
true for the case of molecules, and the strong-coupling 
limit is required to cool efficiently, as we will see in Sec- 
tion IVH 



V. CAVITY COOLING IN THE STRONG 
COUPLING LIMIT WITHOUT SATURATION 



cooling 



cooling hkvY cyV) 



mv^/2 



(22) 



where E is the atom's kinetic energy. From this equation, 
we see that the cooling rate is optimal when the pump de- 
tuning offsets the Doppler shift of a counter propagating 
atom, |Apc| — k^v « 0. When we account for the scatter- 
ing from the counter propagating transverse beam which 
is of detuning |Apc| + k^v, we find the optimal pump- 
cavity detuning to be |Apc| ~ k + k^v. In other words, 
the pump beams should be chirped towards smaller de- 
tunings as the velocity distribution is compressed. 

Like free-space Doppler cooling, the velocity capture 
range of cavity cooling is limited. From Eq. |20] we can 
see that the capture range is limited to ~k about a 
center frequency, kvQ. The maximum cooling force on 
an atom moving in the 2D plane is |27j : fmax ■ vq — 
—HTe{vo)\(kx T kz) • vo|. By balancing the cooling rate 
with recoil heating, one obtains a temperature limit pro- 
portional to Tf (X k(1 -I- C~^). Note that for narrow 



In the good cavity limit ([7_L,.g] > k) or, more re- 
strictively, the strong coupling regime {g > [k,^±\), the 
perturbative treatment described above becomes less ap- 
plicable since an intracavity photon can be reabsorbed 
coherently by the atom many times before finally being 
dissipated via the cavity mirrors or spontaneous emission. 
The atom and cavity system can no longer be treated in- 
dependently and it is advantageous to study the system 
in the dressed-state picture as shown in Fig. [sj In this 
figure, the state space of the first excited states of the 
joint system, |±), is being probed by the pump laser of 
Rabi frequency f2p. The addition of the secondary ground 
state 1 6,0) accounts for all of the Raman loss channels 
found in r ealis tic molecules, and is discussed more fully 
in Section 



VI 



The inclusion the factor of V-/V for the 
N copies of the intracavity particles will be discussed 
in the context of multiparticle effects in Section |VIII| 
The pump laser can only excite the |e, 0) level since the 
transition to |a, 1) is not electric dipole allowed. At low 
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|a,0> = |g> 



FIG. 5: Dressed-state spectroscopy for N three-level atoms 
(molecules) coupled to a cavity mode in the presence of 
a transverse pump field, Qp. The detunings are set to 
[Apa , Ape] < for cavity cooling, and the pump field con- 
nects the absolute ground state of the joint system \g) with 
|e, 0). For efficient cavity cooling, the magnitude of f2p is set so 
that the population of |e, 0) is low, suppressing both incoher- 
ent scattering to \g) and Raman scattering to the secondary 
ground states, collectively represented by |&, 0). With these 
detunings, the atom(s)-cavity coupling, g, mixes the excited 
states to produce the dress-states |±). The lower dressed- 
state |— ), may be viewed as a polariton excitation, as it is 
mostly comprised of the cavity excitation. 



saturations, s < 1, this state decays via coherent elas- 
tic Rayleigh scattering, 'jRy, to the ground state, which 
leads to heating. Additionally, inelastic Raman scatter- 
ing, jRn, may occur, which depopulates \g). The coher- 
ent atom-cavity coupling exchanges population from the 
atomic excited state |e,0), to the cavity's excited state 
|a, 1), at the "Rabi" frequency 2g. Thus, population can 
be decoupled from atomic decay if 5 > 7_l, and \a, 1) will 
decay via mirror leakage at the rate 2k, which leads to 
cooling. 

An alternative, but wholly equivalent method for de- 
riving Eqs. [16] and [17] is to view the pump field as spec- 
troscopically probing the dresses states |±), which are 
found by diagonalizing the Hamiltonian in Eq. [T] Under 
the conditions of Eq. [15] the states are: 

1+) = ce|e,0) +Ca|a,l) (23) 

I-) = Ca|e,0)-Ce|a,l), (24) 

with 

|cep = l-g'/Al,, (25) 

|c,p = gVAl^. (26) 

The decay and populations of |±) may easily be 
obtained — with the justifiable assumption that the coher- 



ence between them vanishes rapidly — by treating them as 
independently pumped by flp. 

Energetically, the cooling arises in a manner similar 
to that described in Section |IV| but whereas the sign 
of Apa was not crucial to cooling in the weak coupling 
regime |19[ I27j . in the good cavity limit, cooling with- 
out atomic saturation is best when A^^ < as shown in 
Figs. [2] and [5] Furthermore, the damping of motion along 
the cavity axis in this regime is instead due to a cavity 
Sisyphus-like effect which arises from the time-delayed 
cavity field with respect to the atomic motion [55] • Es- 
sentially, the dressed-state potential |— ) is steeper as the 
atom climbs than when falling down the potential, which 
provides damping along the cavity axis. This effect is 
also found when cooling via driving the cavity mode it- 
self [25] [221 131] ; though the role of the detunings is re- 
versed and the scaling of the cooling rate with N is less 
favorable [35] . Additionally, a Doppler-like cooling mech- 
anism along the pump field axis has been identified — 
referred to as "polariton" cooling because the cooling 
arises from driving the |— ) state — in which a strongly- 
coupled atom is cooled when its motion is perpendicular 
to the cavity axis [M]- Together with the cavity Sisyphus 
cooling, damping of the particle's motion in 3D is again 
possible as in the perturbative case. As noted in Ref. [37] . 
the capture range can be enhanced in this regime by a 
factor of nearly 10 if one is willing to have the tempera- 
ture limited by h'-f/ks rather than hn/kB- This is a small 
price as the temperature of ground state, polar molecular 
samples is currently limited to T > IQQh^ /ks- 

Refs. [2B] and [10] review these effects and connect how 
the cooling mechanism evolves from weak to the strong 
coupling regimes. Tractable diffusion and damping ex- 
pressions are obtained in these works by keeping only 
lowest-order terms in kv/n. More recent work extends 
this to all orders in velocity, completing the analysis of 
all the single-atom cavity cooling regimes ^37j ^38, .39^. 

A more sophisticated treatment of the cooling rate con- 
firms the role of cooperativity as a benchmark to whether 
the particle can cool before spontaneously scattering. 
Specifically, the ratio of the average velocity damping 
rate, /3, to the free-space scattering rate is [40] : 

/? ^oj,,, Im{{D*f{zl-g^)] 

where Zc = ~K+iApc, Za = —^+iApa, and D = ZcZa+g'^- 
Figure [6] plots this ratio versus the cooperativity and the 
detuning from atomic resonance. The equations for the 
ratio of cooling to spontaneous emission in the cavity 
pumping case (i.e., fid ^ and flp = 0) are cumbersome 
and not listed here. In both cases, the ratio is less than 
unity when C < 1, as expected. This holds true for 
cooling to the blue of the cavity resonance as well (see 
Ref. [20j for information on blue-detuned cooling). 
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FIG. 6: Ratio of the cooling rate to the spontaneous emission 
rate, /3/(7{(t+(t_)), in which a negative- valued ratio indicates 
cooling. Parameters chosen for an OH molecule electronically 
excited in the cavity described near the end of Section |VII| 
which has a finesse of 5000 and length of 2 cm. The detuning 
Ape is set such that the pump field is — k detuned from the 
lower dressed state. 



VI. THREE-LEVEL ATOMS REPRESENTING 
MOLECULES 

As discussed in the last section, single two-level atom 
cavity cooling is well understood. However, the efficacy of 
these cooling mechanisms comes into question when the 
scatterer has a more complicated internal structure. It 
is the central point of this Article to extend the previous 
sections' analyses to the case of cooling molecules, and 
show that C > 1 is an additional, necessary requirement 
for efficient cooling. 

We may investigate multi-level scatters, such as 
molecules, by bundling all higher states into a single state 
denoted and assume that the pump laser is so far 
detuned from cub that this level can never be depopu- 
lated once a single Raman scattering event has occured. 
In other words, no coherence is developed between |6, 0) 
and either |a, 0) or the dressed-states |±) (see Fig. |5|. 
With this approximation, we need only to add an extra 
decay term to the master equation, Eq. [3] to incorporate 
the third level and, consequently, the set of non-ground 
states from which population effectively never returns: 

p = -i H,p + 7_R,y(2(T_po-+ - o-+(7_p - /9(T+o-_)/2 
"7i?,n [o'+o'-j p] /2 + n{2apa} — a}ap — pa) a). 

(28) 

For atoms such as cesium, only one additional laser is 
required to repump the atom back to its ground state 
where the cooling laser, associated with the pump field 
rip, can continue the cooling process. Molecules, on the 



other hand, need more than one repumping laser to pre- 
vent population loss among the myriad vibrational, ro- 
tational, and spin-orbit levels. Unfortunately, the gen- 
erally low free-space branching ratio, T = ^Ry/^Rn, of 
Rayleigh-to-Raman scattering for molecules results in 
population shelving after only a few photon scattering 
events, thereby prematurely quenching the cooling pro- 
cess. This ratio is fixed by the molecular structure, and 
it cannot be modified by detuning except when using 
detunings that are incredibly large Apa « uj^ |31j . In 
other words, even though the molecule is unsaturated 
and coherently scatters the pump field, there remains a 
fixed — and generally large — probability for Raman loss 
per Rayleigh scattering event, regardless of how large 
the detuning. Appendix |B] lists the candidate cooling 
transitions for the radical OH, along with the respective 
T's. 

While building several repumpers becomes pro- 
hibitively complex and expensive, one may use a very far- 
off resonance transverse pump laser to simultaneously ad- 
dress all higher lying metastable states, as was described 
recently in Ref. 41J. However, doing so naturally de 



creases the coohng rate by the ratio fip/Ap^ as in Eq. 16 
This low cooling rate in the far detuned case makes prac- 
tical implementation difficult because of the short inter- 
action time between hot molecules and the small cavity 
mode volume. One can attempt to compensate by us- 
ing a multimode cavity, as explored in Section |VII| and 
Refs. 1^ HI]. The only solution for preventing popu- 
lation shelving is to ensure that there is a vanishingly 
small probability that the molecule will Raman scatter 
during the cooling time. Because T is typically no larger 
than unity for molecules, there is a roughly 50% chance 
that the molecule will be shelved in a dark state after 
each spontaneous scattering event. It is therefore imper- 
ative that the excited state remain sufficiently unsatu- 
rated during the cooling process so that no spontaneous 
free-space scattering is likely to occur. Cavities useful for 
cooling large samples occupy too small a solid angle for 
appreciable Raman suppression. We will discuss these 
practical issues in more detail in Section [Xj 

While T cannot be modified (even by large detuning), 
the rate of scattering into the cavity, and thus the cool- 
ing rate, can be made to be larger than the Raman scat- 
tering rate by a factor equal to (1 -I- T)C. The ratio 
T is for most molecules no larger than approximately 
unity. Therefore, the rate of scattering into the cavity 
versus Raman loss is approximately equal to the coop- 
erativity, C. In other words, with respect to the lower 
dressed-state in Fig. |5] the cavity decay channel (cx 2k) 
can be made dominant over the atomic decay channels 
(c< {7_Ri/, 7_Rn}) if only one can achieve C > 1. Raman 
scattering is not simultaneously enhanced by the cavity 
because its frequency is far-detuned from the cavity res- 
onance and does not spectrally fit into the cavity band- 
width. This situation is markedly different from that 
of cavity cooling two- level atoms, where efficient cooling 
could be achieved even if the cooperativity is less than 
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unity. The cooling rate is slower, but the atoms are not 
shelved into dark states as occurs with molecules. This is 
the essential difference between cavity cooling atoms and 
molecules, namely, molecules require C ^ 1. In the fol- 
lowing sections, we will discuss methods to enhance the 
cooperativity for practical implementations of molecular 
cavity- assisted laser cooling. 



VII. ENHANCING SINGLE MOLECULE 
COOPERATIVITIES WITH MULTIMODE 
CAVITIES 

One method that can be used to achieve larger coop- 
erativities is to decrease the volume of the cavity while 
simultaneously increasing its finesse. This is the route 
taken by single atom cavity QED research [29, 30 , but 
is not suitable here as one wants to cool large diffuse 
samples several millimeters in width and would not fit 
into the required sub- millimeter sized cavities. More- 
over, to maintain a high /k ratio, one needs to increase 
the finesse to the 10'* or 10^ regime. Unfortunately, the 
strongest molecular electric dipole transitions (in light 
molecules like OH) are typically in the blue to UV wave- 
lengths, and cavity finesses much greater than 10^ are 
currently commercially unavailable in the near UV. 

The geometry of the cavity is primarily constrained 
by the molecular beam size or trapped cloud dimension 
and by the inability to obtain high-finesse mirrors in the 
UV. Moreover, there is a trade-off between cooling vol- 
ume and g: a large cavity waist provides a large cooling 
volume, but g will be correspondingly smaller. In the 
case of OH produced by a Stark decelerator, the cavity 
length must be at least > 5 mm to ensure that most of 
the molecules are enveloped by the mirrors, but < 20 cm 
to fit inside typical vacuum chambers. For OH excited on 
the Pi(l) transition (see Appendix|B]) , C = 9.4x 10^^ for 
the TEMqo mode of a cavity of radius of curvature (R) 
and length (L) approximately equal to 2 cm and finesse 
F = 5000. This cooperativity is certainly not sufficient 
for ensuring that the cooling rate will dominate the Ra- 
man loss rate. 

One may wonder whether cooling on purely vibrational 
transitions is more feasible given the high finesse coatings 
available in the infrared and the relatively few number of 
repumping lasers required to close the transition. Unfor- 
tunately, the low decay rate of the vibrational transitions 
mitigates their utility for realistic cooling in all but long 
lifetime traps (see Section|X|. For instance, OH's first vi- 
brational transition at A = 2.8 /im, which naturally has 
the largest T — equal to 1.6 — has a slow, j — 2% ■ 2.7 Hz, 
decay rate |42]. With a single mode R L = 1 cm, 
F = 10^ cavity, the cooperativity is C = 34, and the 
best achievable rate of scattering into the cavity would 
be C • 7 = 27r • 100 Hz. The first vibrational overtone, 
at 1.4 ^m, possesses a lifetime roughly a factor of two 
smaller, but has a wavelength at which it is much easier 
to obtain high power lasers. A similar cavity would give 



C = 17 and a maximum scattering rate into the cavity 
of 2tt -23 Hz. However, in both cases the cavity waist 
is only ~60 /im, which decreases the time fast molecules 
would spend in the cavity mode. 

Before examining how multimodc cavities can help in- 
crease C for electronic transitions, we note that increas- 
ing the single-particle cooperativity by "seeding," i.e., 
driving the cavity with a nonzero fl^, does not help to 
stimulate more photon exchange from the pump beam 
to the cavity mode. The same number of photons would 
be stimulated back from the cavity to the pump beam, 
thereby canceling the energy loss and adding recoil heat- 
ing and atomic saturation. 

Given these restrictions on minimum cavity length and 
maximum finesse, another method for increasing single 
molecule cooperativity is to increase the number of cav- 
ity modes available. One could do this by wrapping more 
cavities around the molecule, but this is highly imprac- 
tical. Near-degenerate cavities — such as confocal, near- 
planar, and concentric — offer an ideal solution |27j in that 
they can support many modes within the cavity linewidth 
K. The molecule can now scatter blue-shifted photons 
into many modes, effectively increasing 5 by a large mul- 
tiplicative factor, rieff . 

Before estimating this factor, we first incorporate a 
multimode cavity into the master equation: 



P = 



- 7-Rn [o-+o--,p] /2 



M 



H 



' ^ Ki(2ai/9di'^ - a^^dip - pa^^ai), (29) 

4 = 1 

M 

a+a- — Ape '^ala, + np[a^+&+]/2 



4 = 1 



M 



4 = 1 



d|o-_ -I- a+Qi 



rf[a. + al]/2}.(30) 



We assume that all M cavity modes within the band- 
width K couple to the atom so that we may remove the i 
indices — effectively treating the modes as a single super- 
mode — and take X^t^i 9i — '^eff 9- However, the coupling 
may not be equal for each mode at each position of the 
scatterer. We next examine how the magnitude of n^g 
may be calculated. 

Vuletic et al. |27j and others previously have shown 
that confocal cavities provide a better compromise be- 
tween cooling rate and cooling volume over the other 
cavities (though the concentric has superior cooling rate 
performance |27j). Consequently, we will focus solely 
on the confocal cavity geometry. For this geometry, 
K — t:c/2RF, where the length and radius of curva- 
ture are equal, L — R, and the cavity waist is simply 
wq = \fRjk j44j. Realistic cavity mirrors have spherical 
aberration which limits the gains otherwise achievable 
with a confocal cavity. Nevertheless, as long as F < fci?, 
using a confocal cavity does increase the cooperativity 
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over the single mode case. In the UV, with the largest 
finesse one could hope to obtain, F w 5000, the cavity 
length must be less than 1 mm for a single-mode cavity 
to be optimal, which is too small for accommodating the 
diffuse molecular cloud. 

The cooperativity of a single-mode cavity, C — 
2FAri/7r = GF/irkL, can be enhanced to the following 
value, which is limited by spherical aberration (sa) (as- 
suming a dipole scatterer oriented perpendicular to the 
cavity axis): 



Csa = 3FAr!,„/4^2 ^ ^^2F/7TkR. 



(31) 



The following ratio provides an estimate of the effective 
number of additional modes: 



"eflf 



~c~ 



nkR 
~2F 



1/4 



(32) 



In the above, R is equal to the length of the confocal 
cavity and F is the cavity finesse. The solid angle of the 
confocal cavity — which is still much less than unity for 
cavities of interest — may be related to that of a single 
mode cavity of equal length and F by: 



Ar2,„ = 8neffAr2/3. 



(33) 



For F = 5000, i? = 10 cm, and A = 308 nm — parameters 
suitable for first electronic transition in OH — the effec- 
tive enhancement is rigff = 5. As a concise figure of merit, 
one would like to have a confocal cavity that simultane- 
ously maximizes spherical aberration-limited cooperativ 
ity and cooling area near the cavity waist, A « ttw, 
product is 



l The 



Csa-A — 



1%FR 



(34) 



We see from this expression that the optimal mirror 
quality and geometry is sensitive only to the product of 
F and R, and a longer confocal cavity is favorable for 
fixed finesse. This expression is not entirely fair, how- 
ever, because a confocal cavity's mode volume is not lim- 
ited by the TEMqq waist, wq, but rather the convolved 
waists of all the accessible modes. For a spherically 
aberrated confocal cavity, this waist is Wsa — '^Rsa = 
2{2TrR^ /kF)i [27 , whose ratio to Wq is: 



Wsa 
Wo 



2{kR/F)' 



(35) 



This is always greater than unity for useful confocal cav- 
ities. 

In addition, these figures of merit do not factor in cool- 
ing time, proportional to k (see Eq. 16 1, which is im- 
portant for a non-stationary molecular sample. We will 
discuss this consideration more in Section |X] Wavelength 
is not generally a tunable parameter, but redder transi- 
tions are favorable. For the Pi(l) OH transition listed 
in Appendix |B] a L = i? = 2 cm, F = 5000 cavity gives 



[Csa, K, 90, Wo, Wsa] = [1-1, Ztt • 7.5 X 10^ Hz,27r • 9.0 x 
10^ Hz, 30 fj,m, 0.3 mm]. Using the confocal cavity, a fac- 
tor of ^12 has been gained versus the single mode case 
(for which C = 0.09), but even with the aid of a confo- 
cal cavity, one still cannot achieve a cooperativity much 
greater than unity for cavities accommodating samples 
of OH. For a, R = L = 10 cm cavity, which would be bet- 
ter for molecular sample insertion and which also may be 
easier to obtain at a finesse as high as F = 5000 in the 
UV, the cavity parameters become [Csa, go,wo, Wsa] — 
[0.47, 27r • 1.5 x 10^ Hz, 2tt ■ 1.8 x 10^ Hz, 70 /im, 1 mm]. 

In practice, this Csa is most likely an upper bound. In 
the experiments of Vuletic et al. jlSJ HHl HZ]) the cavity 
happened to be misaligned from perfect confocality by 
~20 fim in the plane perpendicular to the cavity axis. 
This splits the degeneracy of the cavity modes over a 
bandwidth of 200 MHz from the TEMqo mode position. 
Consequently, for their R — 7.5 cm, F = 2000 near- 
confocal cavity at 852 nm, ji^g should have equaled ~4.5, 
but was experimentally found to be between 2.7 to 3.2, 
or 60-70% of the expected value at the frequency of max- 
imum mode density [47J . If we take this reduction as a 
pessimistic bound, then one can expect a cooperativity 
of Csa ~ [0.73,0.31] for the OH cavity cooling system of 
cavity length 2 cm and 10 cm, respectively. The mode 
volume was 200 times larger than that of the TEMoo 
mode [47j , which is 95% of what Eq. |35] would predict for 
an ideal, spherically-aberrated confocal cavity. 

The multimode enhancement, neff, may be more accu- 
rately calculated by taking into account the actual con- 
focal mode structure, as noted in Ref. [55]. This is done 
analytically using modes that are effectively uniform over 
the cavity length and by employing a numerical calcula- 
tion that includes the Gouy-phase term. This results in 
the following expression for rigg: 



(2Af' + 1)! 
{2M')\\ 



(36) 



where 2M' is the maximum mode index and M — 
(M' + 1)^ are the total number of modes seen by the 
scatterer. Until the paraxial approximation breaks down 
for very large M, i.e., when M\ ^ L, rigg grows roughly 
linearly with mode number as we expect from the pre- 
vious discussion. The number of modes supported by 
the cavity of Vuletic et al. [17] was measured to be 220, 



which produces rigff = 4.3 from using Eq. 36 This is con- 



sistent with Eq. 32 s value of 4.5, and the experimental 
near-confocal cavity realizes ~70% of this enhancement. 

For predictions of future cavity performance, Eq. [32] 
seems to be sufficient. Using the numerical confocal cav- 
ity calculations, the authors of Ref. [3S| reveal a construc- 
tive intracavity mode interference effect that reduces the 
temperature limit, hn/kB, by as much as 20% for scat- 
terers offset from the cavity midpoint. The cooling rate 
could be increased as well, though this effect would need 
to be confirmed with future simulations. While multi- 
mode cavities can aid in increasing cooperativity, exper- 
imentally realistic confocal cavities cannot push C much 
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beyond unity. 



VIII. COMPLICATIONS AND BENEFITS OF 
MULTI-INTRACAVITY SCATTERERS 



As we have seen in the previous section, efficient molec- 
ular cavity cooling is not possible when C < 1, as is the 
case for the Pi(l) transition of OH when the cavity is 
made long enough to accommodate an experimentally re- 
alizable sample. The Raman loss rate will dominate the 
cooling rate, prematurely quenching the cooling process. 
For multiple intracavity molecules, cooling is still impos- 
sible as long as the molecules act independently of one 
another and is unmodified by many-particle effects. In 
other words, because the Rayleigh-to-Raman scattering 
ratio, T, will never be much larger than unity with no 
repumpers (or even with a single one) , the C > I regime 
is necessary to enhance coherent Rayleigh scattering into 
the cavity mode over the Raman free-space scattering. 
To achieve this strong coupling in the absence of collec- 
tive effects, one has no other recourse but to make small 
cavities, an impractical compromise if one intends to cool 
large samples. 

The situation of multiple intracavity scatterers adds 
many-body complexity to the cooling physics, but also a 
means to increase the pcr-molccule cooperativity beyond 
unity. To understand how this collective enhancement 
may arise, let us first consider TV randomly spaced scat- 
terers (which we will presently refer to as atoms) inside 
a cavity, as depicted in Fig. [Tj For simplicity, we assume 
all TV atoms are located at the antinodes. Atoms at the 
nodes do not couple to the field {^p — 0) and those offset 
from both the node and antinodes have a diminished cou- 
pling (■0 < 1) . We can renormalize N to N^.^ = V'l 
and consider the N^^g atoms equally distributed among 
the antinodes (we will drop the eff subscript from this 
point forward). There are two antinodes per wavelength 
and the electric field oscillates exactly tt out of phase 
between them. We designate every other antinode as 
"even" or "odd" to highlight this phase difference. Fig- 
ure |7]depicts the even and odd antinodes, along with the 
on average 7V/2 atoms in the even sites and 7V/2 in the 
odd. In the absence of motional fluctuations, there are 
exactly the same number of atoms at the two types of 
antinodes. For every "even" atom that scatters a photon 
in the cavity mode from the pump beam, there will be 
another "odd" atom that scatters a A/ 2-displaced photon 
along the cavity axis. These two photons destructively 
interfere, preventing cavity field build-up. 

The atoms have a finite temperature, however, and sta- 
tistical fluctuations will cause momentary imbalances in 
the particle number at the even versus odd sites. Assum- 
ing Poissonian fluctuations about the mean, at any given 
moment there are {N + \/N)/2 atoms in the "odd" sites 
versus (iV — ^/N) /2 in the even. Of course, the choice of 
excess atoms in the odd sites is arbitrary, but we will as- 
sume this for concreteness. These additional \/N atoms 
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FIG. 7: a) Cartoon of multiple intracavity particle scatter- 
ing and self-localization. Green dots represent the scatterers 
and the blue sinusoidal line is the intracavity field, E. There 
are two antinodes per wavelength, distinguished by opposite, 
even or odd phases. The even are represented as the upper 
antinodes while the odd are the lower, b) Below a thresh- 
old transverse pump Rabi frequency, f2p < Q.th, the atoms at 
the even and odd sites experience an identical dipole trapping 
force, proportional to \E\'^ . The random nature of the atom 
position implies that at any given time there could be \/iV 
more atoms in the odd sites than the even, c) This leads to 
an increase in scattered light from the unpaired "odd" atoms, 
and the interference with the pump beam creates an ever 
deeper optical potential for atoms in the odd sites, d) The 
even sites are depopulated over the odd, leading to a phase 
transition from particles located every antinode to a A spac- 
ing. Concomitantly, superradiance ensues, which increases 
the per atom cooperativity by TV. The choice of odd over 
even is arbitrary and the symmetry is spontaneously broken 
in favor of one versus the other. 



in the odd sites are unpaired by any atoms located at in- 
teger multiples of A/2 and therefore scatter photons into 
the cavity without destructive interference. Thus, a cav- 
ity field is built-up with Rabi frequency proportional to 
\/N g. This produces a per atom cooperativity equal to: 



1 i^Ng) 



N 



K7 



51 
K7 ' 



(37) 



which is equal to the single particle cooperativity. 

No multiparticle cooperative effect is seen in the cavity 
scattering rate — which is A^F^ — even though the spec- 
trum of the joint atoms-cavity system exhibits a eigen- 
mode splitting proportional to VNg. This highlights the 
difference between pumping the cavity mode itself (fid) 
versus pumping the atomic mode directly with a beam 
transverse to the cavity axis (f^p). In the latter case. 
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the atoms act as independent scatterers when observ- 
ing the rate at which photons are coupled to and then 
leaked from the cavity mode. In the former, the atoms 
collectively act as a giant dipole that modifies the cavity 
transfer function (atoms-cavity spectrum). 

The situation is dramatically modified if the pump is 
made stronger than a critical field Hp > ^Ith [Ml l49] . 
A more detailed explanation will follow after this gen- 
eral overview. Because the cavity field is red-detuned 
from the atomic resonance, an optical dipole trap may 
be formed with trap minima centered at the cavity antin- 
odes. However, an intracavity field will not be formed if 
there is no population imbalance between even and odd 
sites. As the particle positions fluctuate, a population 
imbalance will form and a cavity field will be generated. 
If this thermal position fiuctuation is weaker than the 
ensuing optical dipole trap, then the imbalance between 
even and odd wells will grow in a runaway (positive feed- 
back) process: Eventually all the atoms will migrate to 
favored set of wells as this further increases the cavity 
field and the trap depth. A lattice of periodicity A will 
form in every other antinode of the intracavity dipole trap 
whose strength is proportional to \E\'^ (see Fig.^. Since 
scattered pump photons into the cavity mode now con- 
structively interfere, a large cavity field is built-up pro- 
portional to Ng. The resulting collective single-particle 
cooperativity is now 



cal cavity field: 



1 {Ngy 



N 



(38) 



which is N times the single-particle cooperativity, Cn — 
NC. The total cavity emission is now superradiant, 
TV^Fc, [Sni [5T] and the per particle scattering rate into 
the cavity mode is collectively enhanced to NVc- Since 
one can place many molecules inside the cavity mode, the 
collective cooperativity can be much greater than unity. 
Cat ^ 1. Raman scattering can now be completely sup- 
pressed relative to the elastic scattering rate. 



IX. THRESHOLD FOR SUPERRADIANCE 



The question remains how exactly to trigger this phase 
transition from a random particle distribution to one 
of self-organization. Note, this is a spontaneous sym- 
metry breaking transition — initial fiuctuations determine 
whether the particles localize at the even or odd sites. 
We now revisit the formalism of Ref. [49] to better un- 
derstand the problem. From a Hamiltonian that incorpo- 
rates the weak driving approximation, Eq. |15| we obtain 
the stochastic time evolution equation for the semiclassi- 



Apc - C/o^cos(fczj) 

] 

K -I- Fq y~^cos^(fczj) 

3 

77(,f[ cos(fczj) cos{kxj) + ^Q,, (39) 



where Zj and the positions of the jth atom along 

the cavity field and pump field, respectively, and is 
the Langevin noise term for the cavity field. Uq and Fq 
are the dispersive drive and decay rates of the cavity field 
via the atomic medium, respectively: 



Un = 



Fn = 



+7i 



2 • 



(40) 



The effective cavity pump strength due to coherent scat- 
tering of the pump field by the atom is: 



-iA 



pa+l± 



Expressions rj^^/ ^ and Fq reduce, respectively, to Eqs 



(41) 
17 



and 
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when |Apc| ~ k,- 



The force on the atoms along the cavity and pump 
axes, respectively, are: 



Pz 



d 

-hUo\a\^ ^— cos^ikzj) 
dzj 

d 

ih{'r]*ffa — ileSC(*))j) — cos{kxj) cos(fczj) + 
ozj 

d 

—hUo{Hp/g)^— — cos^{kxj) 



dxj 



d 



ih{r]lQa - 77effa*)) ^ cos{kxj) cos{kzj) + , 



(42) 



where and $,x- are Langevin noise terms for the motion 
along the z and x, respectively. If we look at the sign of 
the interference term between the cavity and pump fields, 
proportional to r]gg, a checkerboard-like pattern appears 
in the crossed intracavity and pump fields. A single line 
of antinodes along the cavity axis is depicted in Fig. [7] 
Moving one half-wavelength perpendicular to the cav- 
ity axis — along the pump field and in the x direction — 
changes the phase of the pump field while the cavity field 
phase stays constant: The spacing of the even-odd antin- 
odes are shifted by A/2, thus creating the checkerboard 
pattern. If the even and odd antinodes are equally pop- 
ulated, then the sum over cos{kzj) cos{kxj) is zero and 
there is no contribution to the cavity field from the rygfj 
term. Once the population becomes imbalanced, the cav- 
ity field is pumped by this term and grows in a runaway 
process as the atoms migrate to one set of antinodes. 
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The optical trapping potential felt by each atom along 
a string of antinodes parallel to the cavity axis is 



V{z) = U2 cos^(fcz) + Ui cos{kz), 



(43) 



which may be understood from the first of Eqs. |42| with 
Ui the strength of the interference term between the 
transverse pump and cavity mode and U2 the dipole trap 
from the cavity mode. The depths of the potential are 



U2 = hloUoN^ {cos{kz)Y 



(44) 



Ui = 2hIoN (cos(/cz)) [Acp - NUo {cos^{kz))],{A5) 



where 



lo 



\VcS\' 



[k + NTo (cos2(fcz))]2 + [Acp - NUo (cos2(fcz))]2 ' 



(46) 

This latter expression may be understood as the per atom 
scattering rate into the cavity mode and is a generaliza- 
tion of Eq. 16 for > 1. To understand why one set of 



antinodes are preferred over another, we look at the signs 
of Ui and 1/2- When Uq < 0, as is the case with negative 
detuning from the atomic line, the U2 term represents a 
cavity field that functions as a red-detuned optical lat- 
tice with equal trap depth in the even and odd antinodes. 
It is the Ui interference term that breaks the symmetry 
between the antinodes and gives rise to the checkerboard 
pattern of atoms occupying only every other antinode. 
We describe this process as follows: The spatial average 
of atomic positions along the cavity axis, {cos{kz)), is 
either closer to -1-1 or to -1 initially. If (cos(A:z)) < ini- 
tially, then there are more atoms near the odd sites and 
consequently, the odd sites have a deeper potential than 
the even. More atoms localize in the stronger wells of 
the odd sites, further reinforcing the asymmetry between 
the odd and even antinodes. This results in a runaway 
self-localization process, leading to (cos(fcz)) assuming a 
value closer to -1. The same holds true if (cos(fcz)) is 
initially -1-1, except the atoms eventually localize around 
the even sites. It has been shown that the final state — all 
atoms at the even checkerboard sites or all at the odd — is 
stable once formed [HI HI] . 

We now address the criteria for reaching threshold, 
rip''. The threshold condition has been derived in two 
ways [15] . The first method employs a mean-field approx- 
imation that assumes an intracavity gas of constant den- 
sity, which is valid in the thermodynamic limit N —^ 00, 
g ^ 0, K — const, and Ng'^ cx N/V — const. For driving 
the lower dressed state below resonance by an amount 
equal to K, i.e.. Ape = -^C^o ~ where NUq is the energy 
shift of the I—) dressed state, the threshold is: 



^7) 



knT k\A 



pa 



Ng 



■V2. 



(47) 



We see that the depth of the optical dipole trap at thresh- 
old scales linearly with temperature {Vth oc f^^^), is scaled 
inversely by the temperature limit Tf « hn/kB, and is 



inversely proportional to the rate at which pump photons 
are scattered into the cavity mode, cx Ng"^ /kA^^. This 
may be simply understood from the statement that the 
trap depth must be large enough to quench the diffusion 
due to thermal energy. 

A threshold that scales inversely with ^/N would not 
necessarily prevent the triggering of self-localization in 
experimentally realizable samples of molecules, but nu- 
merical calculations [H] indicate that a hysteresis effect 
in the phase transition pushes the onset out to larger 
pump fields. The numerically verified threshold becomes: 



f2p > flth — 



ksT K\Apa\ ^TT 

hn m/^g 2 ' 



(48) 



This significantly worse scaling with N would prevent the 
triggering of the superradiance. 

A major impediment to triggering threshold is the need 
to do so without saturating the transition. In other 
words, one cannot simply increase fip without a cost in 
the number of spontaneous Raman emissions. While fip 
must be greater than threshold: 



ksT n\Apa\ 



(49) 



where x might equal to 2 or 4, the saturation condition 
must additionally be satisfied: 



^ <i. 



4A2 

pa 



(50) 



This implies that the molecule number must be much 
greater than: 



kRT 



hn 2g^/s 



(51) 



Using the following parameters for OH in the 2 cm con- 



focal cavity described in Section VII and assuming T— 10 
mK, we have Nq > 8.5 x 10^ for x = 2 and Ao > 7.3 x 10^ 
for X = A. As discussed in Section |X] A'o could be achiev- 
able for X = 2 with improvements in Stark deceleration, 
but not likely achievable for x = A. Thus, the question 
of whether the threshold scales as N~^/^ or as N~^/^ is 
of paramount importance. 

Seeding the cavity by driving the cavity mode itself 
with rid 7^ is one possible method to increase the effec- 
tiveness of the transverse cavity cooling scheme, and the 
phase of the drive field with respect to the pump field has 
been previously shown to affect the symmetry breaking 
of the self- localization process [ID]. Seeding introduces 
an intracavity optical dipole trap without having to first 
scatter from the intracavity atomic medium. One might 
expect this to hasten the cooling process or decrease the 
number of molecules required for threshold, but neither 
effect has thus far been seen in quantum Monte Carlo 
simulations. This may be due to the fact that a lat- 
tice formed by driving the cavity does not contribute 
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to the positive feedback mechanism responsible for self- 
localization. Nevertheless, any possible method for re- 
ducing the threshold requirements for localization war- 
rants more investigation via additional numerical simu- 
lations. 

Reference [55] identified a mechanism that suppresses 
free-space scattering well over what is expected from the 
Cn factor, once self-localization has been triggered. This 
is due to pump-cavity mode interference in which the col- 
lective atomic dipole oscillates out of phase with respect 
to the pump, suppressing the atoms' excitement. This 
effect points to another method for achieving superra- 
diance in the case where N < Nq: increase the satura- 
tion, s, as the molecules enter the cavity to lower the 
criteria for reaching threshold. Once self-localization is 
established, spontaneous emission should be quenched, 
and the process only sacrifices a fraction of the initial 
molecules to achieve self-localization for the remaining 
ones. For example, this technique may be employed for 
cooling on the vibrational transitions which would nom- 
inally require Nq > 10^ for x — 2, but can only Raman 
scatter at a maximum rate 7/2 « 27r x 50 Hz. 

The simulations in Ref. [IS] suggest that the thresh- 
old should scale with particle number as a; = 4, but the 
experimental results on self-localization [46l |47l |53j are 
consistent with x = 2. References [ITJ [53] attribute this 
experimental x = 2 scaling to the onset of recoil induced 
resonances (RIR) [51], but it is noted in Ref. [H] that 
RIR would lead only to transient enhancements as op- 
posed to the long-lived collective states described by self- 
localization. Future work will elucidate the dynamics of 
this phase transition with additional simulations. Ques- 
tions to be explored are: how the transition threshold 
scales with intracavity particle number; the role of seed- 
ing; the characteristic onset time for superradiance; and 
how many molecules are lost before spontaneous emis- 
sion is suppressed in the self-localized state. An addi- 
tional scenario for simulations to explore is the case in 
which the recoil energy of the particle of interest, hujrec, 
is comparable to the final cooling temperature, hn/kB- 
All these questions are of crucial interest for understand- 
ing the efficacy of using transversly-pumped cavities to 
cool molecules. This latter question is of particular in- 
terest for the relatively light OH molecule whose recoil 
energy, 12 /iK, is comparable to hn/kg — 7 jiK for the 
10 cm confocal cavity described above. This naturally 
leads us to the last main section of this Article, exper- 
imental considerations for the cooling of molecules, and 
in particular, OH. 



X. EXPERIMENTAL CONSIDERATIONS 

Most of the physics to be studied with polar molecules 
requires them to be in their rovibronic ground state. 
Supersonic expansions of the molecules entrained in a 
buffer gas readily produces packets of polar molecules 
with quenched rotational and vibrational motion. The 
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FIG. 8: Approximate experimental OH Stark-decelerator ef- 
ficiency curve — density versus final velocity — for the JILA ap- 
paratus [141 1551 156) . Inset shows the approximate volume of 
the OH packet in zones II and HI. The packet is larger in the 
pre-Stark decelerator zone I, and expands to fill the magnetic 
or electric trap volume in zone IV. 



technique of Stark deceleration [ll, can readily produce 
slow packets of these ground-state molecules, and in par- 
ticular, the polar molecule OH [12,. OH is produced 
either via water discharge or photolysis of nitric acid. 
When entrained in Xe or Kr, nearly all of the OH — 
upon expansion — is in the lowest rovibrational ^H3/2 
A-doublet ground states. Although the packet that is 
formed has a much narrower velocity spread than ex- 
pected from a Maxwell-Boltzman distribution, the mean 
center-of-mass (CM) velocity is approximately 400 m/s 
when using Xe as the buffer gas. The packet has a 
^^15% and a ~7% velocity spread in the longitudinal and 
transverse dimensions, respectively. Thus, in the mov- 
ing frame, the molecular packet is cold (on the order of 
1 K). Before slowing the packet to near-zero CM veloc- 
ity, the packet must pass through a skimmer to prevent 
backscattering from collapsing the supersonic expansion 
and to limit the flux of unwanted gas into the decelerator 
chamber. 

A major practical difficulty with cavity-assisted laser 
cooling of molecules involves ensuring that a sufficiently 
large molecular sample experiences the mode volume of 
the cavity for a long enough period of time to cool. At 
best, cavity waists are no larger than a millimeter, and 
ground state polar molecules have yet to be produced 
in sufficient number at CM velocities below '^lO m/s. 
For the technique of Stark deceleration. Figure [8] shows 
the inherent trade-off between slowed packet density and 
final velocity. 

The efficiency curve is divided into four zones, each 
zone being amenable to a different regime of cavity cool- 
ing. Post-skimmer, the OH packet is typicallyof density 
10* to 10^ cm"*^, as shown in zone I in Fig. [s] Because 
the packet is moving at high velocity, it experiences the 
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FIG. 9: a) Transit time for crossing perpendicular to the 
cavity waist. Black (solid) curve for crossing all the modes of 
a confocal cavity, blue (dashed) for crossing only the TEMqq 
mode of the cavity, b) Time to cross 66% of the cavity length, 
L — 10 cm, parallel to the cavity axis. 



cavity mode for only a short time. However, the number 
of OH in the cavity mode may be sufficiently large to 
trigger rapid cooling. Self-localization happens no faster 
than l/n, but seems to occur within tens of microseconds 
or less [46*, '4F. Figure [9] shows the time required to cross 
the cavity waist for a given velocity. It may be possible to 
transversely cool the fast OH packet with a cavity. Since 
Stark deceleration would be much more efficient with an 
OH packet of narrower transverse velocity spread 56J, 
the cooling in zone I could lead to the Stark slowing of 
many more molecules in zones II-IV. 

Zones II and HI offer a compromise between particle 
number and cavity transit time. In zone II, it may be 
possible to increase the molecule-cavity interaction time 
by redirecting the molecule along the cavity axis with 
an electrostatic guide. The molecules are slow enough 
in zone HI to be stopped and reflected by an electro- 
static ^7] or magnetostatic [T3] mirror. During the re- 
flection time, the molecules could spend more than a mil- 
lisecond in the mode of an optimally situated cavity. 

Magnetic [14] and electrostatic [58] traps and AC elec- 
tric traps [591 have been used to trap polar molecules at 
the terminus of a Stark decelerator, in zone IV. How- 
ever, room temperature blackbody radiation limits the 
OH lifetime in these traps to < 3 s J60]. Thus, the cool- 
ing must occur on a time scale much faster than this. 
While 1 kHz cooling rates, for example, are fast com- 
pared to this lifetime, it must be noted that the fraction 
of the trap volume occupied by the cavity waist can be 
quite small. Consequently, the molecules do not spend 
a large amount of time being cooled, and a much larger 
cooling rate than expected is necessary. 

Future work-using numerical cavity cooling simula- 
tions incorporating molecular motion-will explore opti- 



mal cavity and decelerator geometries for cavity cool- 
ing OH. The ground state polar molecular samples are 
currently 1000 times hotter than the cooling limit im- 
posed by the cavity linewidth. Simulations will address 
the efficiency and scaling of cooling such hot samples as 
well. Future improvements to the Stark decelerator tech- 
nique [56J will produce higher densities of slow molecules, 
seemingly an important step towards successfully achiev- 
ing the cavity-assisted laser cooling of molecules. The 
use of other techniques, such as feedback [61] and optical 
Stark deceleration [55] may be necessary for increasing 
the cooling rate and obtaining a larger number of slowed 
molecules. 

In summary, we have identified important necessary 
conditions for cavity cooling ground state molecules, 
whose open channels and high initial temperatures pose 
unique challenges. We determine that to prevent Ra- 
man loss and thereby achieve efficient cooling, the coop- 
erativity should be greater than unity. Several methods 
for increasing the cooperativity and cooling rate are ex- 
amined. These include the use of multimode cavities; 
simultaneously driving (seeding) the cavity mode while 
transversely pumping the atomic medium; and inducing a 
self-localizatizing phase transition of the molecules' posi- 
tions. While multimode cavities are useful for increasing 
the cavity cooling volume and raising the cooperativity 
to near unity, only by inducing the superradiant phase 
transition can one achieve molecular cavity cooling with 
certainty. We have assessed the feasibility of triggering 
this superradiant state as well as address the possibility 
of cooling under current experimental constraints, such 
as UV mirror coating technology, with particular empha- 
sis given to the present performance of OH Stark de- 
celerators. More in-depth numerical simulations of the 
seeding and self-localization dynamics in the presence of 
molecule motion — either from a beam or due to a har- 
monic trap — are required to fully access the feasibility 
of efficient cavity-assisted laser cooing of molecules. Fu- 
ture work will address these questions via both quantum 
Monte Carlo simulation 20J and experiment. 
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APPENDIX A: ATOM-CAVITY COUPLING 
PARAMETERS 

This section describes the various parameters encoun- 
tered when analyzing a cavity QED system and attempts 
to connect this language with that used by other au- 
thors. Single atom strong coupling requires the atom- 
cavity coupling, go, to be much larger than both the 
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TABLE I: Comparison of the three best OH electronic cooling transitions. Pi(l) has a better bare Rayleigh-to- Raman ratio, 
T, but requires one microwave pumping stage and does not have a cycling hyperfine transition. The Qi(l) transition has a 
smaller T and a cycling transition on F" = 1 — > _F' = 2 if two microwave pumping stages are used to prepare the OH in the 
F" — 1 hyperfine ground state. A compromise transition is Q2i(l), which has an intermediate T but no cycling hyperfine 
transition. Note that the Pi(l) transition is 92% closed with just one repumper on the Pi2(l) line. Transition wavelengths and 
lifetimes are from the software packages HITRAN and Lifbase |64j . 
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[J',N'] 
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T 


# repumpcrs 
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^ fi-wave pulses 
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308.256 


[1/2, 0] 


2.32 


1.43 


2 


0.4% 


no 


1 


Qi(i) 


307.933 


[3/2, 1] 


2.32 


0.28 


4 


4% 
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2 


Q2l(l) 


307.937 


[1/2, 1] 


2.32 


0.65 


2 


1% 


no 






atomic dipole decay rate, 7j_ — 7/2, and the decay rate 
of the cavity field, k. Specifically, the saturation photon 
number, uiq = 7^/2(7q, and the critical atom number, 
iVo = 27j^k/5q, must both be much less than unity. For 
conversions: 



TTC 

AQ' 



TTC 

2Lk, 



AO 
2L' 



(Al) 



where L is the cavity length, Q is the quality factor, and 
F is the finesse. 

The inverse of the critical atom number is known in 
the optical bistability literature as the cooperativity pa- 
rameter C = 1/Nq |55]- Bistability ensues when the col- 
lective cooperativity NC > 1, where N is the intracavity 
particle number. Other authors [T9 | |2 H l27 l l45 l l46l l53] 
choose to describe the cooperativity not in terms of the 
notation of go, etc., but rather in the more experimen- 
tally physical language of scattering, which is applicable 
in the weak coupling regime. This is often less convenient 
when working with the master equation, and the follow- 
ing describes how to convert between the two descriptions 
of the cooperativity. The collective cooperativity may be 
reframed as the product of the number of intracavity pho- 
ton roundtrips {F/tt) and the fraction of scattered light 
captured by the cavity mode, 2 Ail = 2 • 3/(fc^w§), where 
k = 27r/A and wq is the cavity mode waist: 



C = 2Ar2F/7r = 6F/{7Tk^w^) = 77. 



(A2) 



This Purcell factor [55], defined as = SQX^/in^V is 
none other than the cooperativity, and is typically the 
preferred notation in semiclassical treatments. The fol- 
lowing demonstrates their equivalence: 

2 1 ^ 3cA^7± 2LF ^ 3QA3 

Since k ~ 27r/A and Wq cx Vm/iLn), upon substitution 
we have: 



C- 6^ 
irk^WQ 



(A4) 



Similarly, the collective enhancement, Nij, is simply the 
collective cooperativity, NC. 



APPENDIX B: CANDIDATE OH COOLING 
TRANSITIONS 
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FIG. 10: OH level diagram depicting the Pi(l), Qi(l), and 
Q2i(l) electronic transitions. The first excited level A '^E"'" 
and the two X ^H grounds states are shown. For clarity, only 
the Pi(l) decay channels are shown: Pi2(l), Oi2(2), and the 
v' = —> u" transitions. 



The Stark decelerator provides samples of weak-field 
seeking F" = 1, X^U^/^ OH molecules [Hj. F" = 1 is 
spectroscopic notation for the lower energy state of a dou- 
blet {F= 1, J = N + 1/2, while for F= 2, J = - 1/2) 
and is not the hyperfine quantum number, F, which is 
italicized. The double-prime denotes ground state quan- 
tum numbers, while a single prime refers to the ex- 
cited state. The total angular momentum is equal to 
J — L + S + R, where the three vectors are orbital, spin, 
and rotational angular momentum, respectively. The an- 
gular momentum N is defined as N = R + L. 

The majority of the OH ground state population ar- 
rives at the terminus of the decelerator in the z/" — 0, 
F" = 2, e = / symmetry, and p = + parity state. Mi- 
crowaves can easily transfer this population to the e, — 
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lower A-doublet states [231 123- The total angular mo- 
mentum in the F" = 1, X'^\l^/2 ground state is J" = 3/2. 
The orbital plus rotational angular momentum quantum 
number in this ground state is N" — R" + L" — 1. 

Electric dipole transitions are allowed between transi- 
tions of opposite parity that satisfy A J = 0, ±1, and 
A7V = 0, ±1, ±2, where AJ = J' - J" and AiV = 
N'—N". The strongest electronic cooling transitions that 
originate in the F" — 1, ground state are Pi(l), 

Qi(l), and Q2i(l). The spectroscopic notation is defined 
as [57]: ANp,p„{N"), where AiV = [-2,-1,0,1,2] ^ 
[O, P, Q, R, S]. In this notation, if F' = F", then only one 
subscript index is used. 

Tablejljhsts the properties of the ^113/2 to ^^^^2 i^v = 
band) transitions of interest. To cool on the Pi(l) 
transition, a single microwave pulse is needed to drive 
the population to the e state. The other two electronic 
transitions need either one additional microwave pulse 
(Qi(l)) or no microwaves at all (Q2i(l)) for initial state 



preparation. Only Qi(l) offers the possibility of a AF = 
+ 1 transition, which allows the molecule to easily return 
to the same Zeeman hyperfine ground state after each 
Rayleigh scattered photon. The recoil frequency for these 
transitions is: ujrec = 27r • 1.23 x 10^ Hz. 

With two repumpers, one each on the Pi2(l) and the 
012(2) lines, the Pi(l) transition remains 0.4% open 
due to the possibility to scatter to higher vibrational 
levels (i.e., 1/ ^ 0). Closing the Qi(l) transition to 
4% requires — in order of importance — repumpers on the 
Pi2(2), Qi2(l), Pi(2), Qi2(3) lines. To close the Q2i(l) 
transition to 1% requires repumping the Q2(l) and P2(2) 
lines. Figure 10 sketches the relevant energy levels for 



the electronic transitions, but for clarity only includes 
the decay channels from the Pi(l) transition. The first 
vibrational transition, Q|e, may be completely closed 
with five repumpers and is ~60% closed with just one re- 
pumper on the P|/ line, where the notation is AJ{J")e 
and e = e, /. 
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